Average Angular Velocity 

Hanno Essen 
Department of Mechanics 
Royal Institute of Technology 
S-100 44 Stockholm, Sweden 

1992, December 



Abstract 

This paper addresses the problem of the separation of rota- 
tional and internal motion. It introduces the concept of average 
angular velocity as the moment of inertia weighted average of 
particle angular velocities. It extends and elucidates the concept 
of Jellinek and Li (1989) of separation of the energy of overall 
rotation in an arbitrary (non-linear) A^-particle system. It gen- 
eralizes the so called Koenig's theorem on the two parts of the 
kinetic energy (center of mass plus internal) to three parts: center 
of mass, rotational, plus the remaining internal energy relative to 
an optimally translating and rotating frame. 
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1 Introduction 



The motion of a rigid body is completely characterized by its (center of mass) 
translational velocity and its angular velocity which describes the rotational 
motion. Rotational motion as a phenomenon is, however, not restricted to 
rigid bodies and it is then a kinematic problem to define exactly what the 
rotational motion of the system is. This paper introduces the new concept 
of 'average angular velocity' as the solution this problem and discusses some 
applications briefly. 

The average angular velocity concept is closely analogous to the concept 
of center of mass velocity. For a system of particles the center of mass velocity 
is simply the mass weighted average of the particle velocities. In a similar 
way the average angular velocity is the moment of inertia weighted average 
of the angular velocities of the particle position vectors. 

The separation of rotational motion from the internal degrees of free- 
dom of a system is of interest in a wide variety of applications. Among the 
more obvious are vibration-rotation coupling in polyatomic molecules and 
the understanding of biomechanical phenomena involving falling cats, figure 
skaters, springboard divers etc. The basic theoretical work on the subject 
is Carl Eckart's (1935) whose work has been extended and elucidated by 
many authors, for example Sayvctz (1939) and Essen (1979). Biomechanical 
aspects have been discussed by Frohlich (1979) from a multibody dynamic 
point of view, and by Essen (1981) from the Eckart point of view. 

Considering the maturity of the subject one might think that all basic 
theoretical results are quite old, but a letter on the subject was published 
as late as 1989 by Jellinck and Li (1989). They showed that one can define 
an angular velocity vector and separate out a rotational part of the energy, 
for an arbitrary (nonlinear) system of particles, without the use of the rigid 
reference configuration that is needed in the Eckart construction. 

In this paper I present some new developments based on ideas related 
to those of Jellinek and Li. In particular I introduce the concept of average 
angular velocity as the weighted average of particle angular velocities with 
moments of inertia as weights. Especially for the elementary, but pedagogi- 
cally important, case of fixed axis rotation this leads to simple and concep- 
tually powerful results. The relevance of these results to spinning athletes 
and turning cats are briefiy indicated. 

Later sections of the paper treat the three dimensional case. It is shown 
that the separation of a rotational part of the kinetic energy can be done in 
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a way closely analogous to the well known split of the kinetic energy to an 
overall center of mass translational part plus the kinetic energy of internal 
motion relative to the center of mass system. A second split is thus done, 
now to an overall rotational energy plus a remaining part corresponding to 
motions in a rotating center of mass system, the rotation of which is given 
by the average angular velocity. 

Throughout this paper I consider a system of particles with masses 
rrik and position vectors Tk {k = 1, . . . , N). A distance from the Z-axis is 
denoted by p, and a distance from the origin by r, the azimuthal angle is if 
and the angle to the positive Z-axis is 9. 



2 Average Angular Velocity Around a Fixed 
Axis 

The z-component of the angular momentum, L^, is by definition 

N N 

Lz = Yl '^ki^k X Vfc) • = ^ mkixki/k - Vk^k)- (1) 

k=l k=l 

Here the Z-axis has a fixed direction and is cither fixed in an incrtial system 
or fixed in the center of mass of the system. We now introduce cylindrical 
(polar) coordinates p, <fi, and z. In terms of these we have for the position 
vectors Tk = pk Gpk+^k and for the velocities tk = Pk e^^^ +Pk'Pk G<Pk+^k e^. 
This gives us 

(rfc X Vfe) • = [{pk + Zk e^) x {pk e^^ + pk^Pk e^^ + 4 ej] • = pl^Pk, (2) 
so that 

TV 

Lz{t) = J2mkplit)Mt)- (3) 

k=l 

If we now define the average angular velocity of the system, around the Z- 
axis, by 

a;av(t) = Mt))t = J, , (4) 

}^k=l^kPk[t) 

and the (instantaneous) moment of inertia, with respect to the Z-axis, by 

N 

Jz = Y. ^kpl (5) 
k=l 
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we see that we get 

Lz = Jz^e.'v- (6) 

If all particles have the same angular velocity, v?, then, of course, = ^■ 
This happens, in particular, if the system of particles is rigid and rotates 
around the Z-axis, but also more generally for any type of motion that obeys 
(pk = (p with arbitrary and i^. For these cases one gets the standard result 
that Lz = Jzp- 

We can now apply these results to the z-component of the angular mo- 
mentum principle (L = M) in the form 

Lz = Mz. (7) 

The general result is that 

Jz^^s^v + Jzi^w, = Mz. (8) 

If the angular velocity is well defined we can replace cUav by in this expres- 
sion. If we furthermore assume that the body is rigid so that Jz = constant, 
equation (jH)) reduces to the standard result Jz^p = Mz. 

If there is no external moment (or 'torque') with respect to the Z-axis, 
so that Mz = 0, then the ^-component of the angular momentum vector will 
be conserved, Lz = constant, and, in view of equation (jH)), we find that 

Jz{t) uJavit) = constant. (9) 

This then says that a large moment of inertia implies small average angular 
velocity and vice versa. Here it is not assumed that there is rigidity or even 
a definite angular velocity. It is this more general form of the standard text 
book formula, Jzp = constant, that is actually 'used' by springboard divers 
and figure skaters. 



3 The Cat Landing on its Feet 

For some arbitrary quantity 7^ the averaging of equation (@)) can be written 



A question of interest is to what extent the average angular velocity can be 
understood as the time derivative of the 'average angle', {(p)t, of the system. 
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The subscript t on the averaging bracket is meant as a reminder that the 
weights in the averaging, mkpl{t), are time dependent and this means that 
the time derivative of an average will not be the same as the average of a 
time derivative. 

If we take the time derivative of 



The average angular velocity is thus not simply the time derivative of the 
average angle. This is, of course, essential if a cat, dropped upside-down, is 
to be able to land on its feet. Equation Q shows that if cUav = initially, 
it will remain zero in the absence of external torque. The above equation 
reassures one that the cat, nevertheless, can change its average angle by a 
proper combination of angular and radial motions. 

4 The Average Angle Concept 

The concept of an 'average angle' requires some comment. The value of this 
angle will, of course, depend on the direction of the fixed reference direction 
(the X-axis). It will also depend on whether one thinks of ip as going from — tt 
to TT or if it goes from to 2n, when one assigns initial values to the (pk- That 
is, (/?av depends on whether the necessary 2% jump comes at the negative or at 
the positive X-axis, respectively. For a cylinder the initial average angle will 
be zero in the former case and tt in the latter. The actual value of the average 
angle therefore has little physical meaning; its significance comes from the 
fact that it defines a reference direction in the particle system. Then, when 
the system has moved, it will tell how large, on average, the net turn has 
been. 

The time dependence of the averaging naturally vanishes if the radii, pk, 
are constant. It is interesting to note that it also vanishes in the more general 
case when the time dependencies of the (cylindrical) radii are of the form 



Pax = 



(11) 



a simple calculation shows that 





Pk{t) = f{t) dk, 



(13) 
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where f{t) is some (positive) function of time and the dk are constants. The 
average then becomes 



For this case then, when the cyhndrical radial motion is a common 'scahng', 
the time derivative operator commutes with the operation of taking the av- 
erage. Consequently the average angular velocity will be the time derivative 
of the average angle and similarly for angular acceleration. 



5 Konig's Theorem 

The kinetic energy of an iV-particle system is given by the sum 

1 ^ 

T = -Y^rukh-h- (15) 
^ k=i 

If one introduces the center of mass position vector 

R.^f^, (16) 

and then re- writes the position vectors of the particles as 

Vk^R + rl (17) 

the kinetic energy is seen to fall into two parts. One part corresponds to the 
motion of the center of mass of the system while the remaining part is due 
to the motion of the particles relative to the center of mass system: 

T= imR-R + T'. (18) 

Here m = Y^=\ f^k is the total mass and T', which is given by 

T'-\Y.m^i'k-i'k. (19) 

^ k=\ 

vanishes if the particles do not move relative to a reference frame in which 
the center of mass is at rest. The absence of cross terms is due to the 
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fact that the sums, J2k=i ''^k^'k — ^k=i — 0; over the center of mass 
system position vectors and velocities, vanish. This result is sometimes called 
Konig's theorem (see Synge and Griffith 1970), or the 'law of the two parts 
of the kinetic energy', if it is given any name at all. It is this result that can 
be taken a step further according to JeUinek and Li (1989) in the sense that 
T' can be split into two parts, one corresponding to an overall rotation and 
the rest corresponding to the motion relative to the rotating system. This is 
shown in the next section. 



6 The Average Angular Velocity Vector 

The quantities of this section may be thought of as referring to the center of 
mass system but we will drop the primes of the previous section. Introduce 
spherical coordinates {rk,0k,^k) for particle k and corresponding moving 
basis vectors (e^^,, e^t^, ei^^^,). The position vector, r^, of particle k is then 

rfe^r^e^^. (20) 

and the velocity of the particle is 

rfc = ''fce^, X Ffc. (21) 

Here u)k is the angular velocity of the position vector of particle k. It is given 
by 

ujk = (pk cos 9k - (fk sin 9k bq^ + 9k e^^ ^^k^z + h ■ (22) 
The kinetic energy of the A'^-particle system is now 

1 ^ 

^' = 9 51 ^k{rk er, + u^fc X Tk) ■ {h e^, + x r^) (23) 
^ k=i 

^ N N 

= 9 Z] iT^krl + -Y.^k{<J^kX rk) ■ {oJk X rk) (24) 
^ k=i ^ k=i 

^ N 1 ^ ^ 

^oYl ^krl + 7^Y1 Jk<^k ■ Wfc (25) 

^ fc=l ^ k=l 

where Jk is the contribution of particle k to the (instantaneous) inertia tensor 
J of the system. The matrix components of the inertia tensor Jk in the basis 
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(e^^ , eg^ , J are given by 



Jk 





mfcr| 




(26) 



rukri 



Using this one easily verifies that the sum of the kth terms of the sums of 
formula (j25j) gives the usual expression for the kinetic energy of particle k in 
spherical coordinates: = \'mk[r1 + r|(v9|sin6'fc + 9l)]. 

Below we will need to manipulate sums of terms like those in the above 
expression for the kinetic energy. One must then remember that the position 
dependent basis vectors (e^^. , e^^, , e^,, ) are different for each particle (as indi- 
cated by the index k). In order to proceed we therefore return to a common 
Cartesian basis in the expression ^ Yl,k=i Jk^k ■ <^k- The inertia tensor will 
then have the matrix components 







1 







M 




Jk = rrik 














-( 






\ 








rl) 





and the sum of these gives 



J 




XkXk XkUk XkZk 

VkXk VkUk UkZk 

ZkXk ZkVk ZkZk 



yl + zl -XkVk -XkZk 



(27) 



-UkXk 



-UkZk 



(28) 



-ZkVk xi + yi 



i.e. the usual (instantaneous) inertia tensor for the system of particles. 
We now define the average angular velocity vector, u^av, by 



N 

Jk'^ki 

k=l 



(29) 



as the inertia tensor weighted average of the individual particle angular ve- 
locity vectors. Here it is necessary that J is invertible so that one can solve 
for ct^av by multiplying to the left by J~^ . This means that the particles of 
the system may not lie on a line since then the inertia tensor is singular. 

If we now denote by uj'k the angular velocity vector of particle k relative 
to the reference system that rotates with the average angular velocity we 
have 

^k = t^av + ^'k, (30) 
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since angular velocity vectors are additive. These relative angular velocities 
fulfill 

N 

E Jku^'k = 0, (31) 

fe=i 

so a calculation completely analogous to that leading to Konig's theorem 
gives 

r = i Ju^av ■ u^av + T". (32) 

Here 

1 Af 1 iV 1 TV 

= ^ E m^rl + ^ E Jku^i ■ ^1 = ^ E "^/^v;: . v^' (33) 
^ fc=i ^ fc=i ^ fc=i 

is the kinetic energy relative to a reference frame that rotates with the average 
angular velocity (around the fixed center of mass) and v^' is the velocity of 
particle k as measured in this frame. 



7 Remarks on the Conservation Laws 

One notes that formula ()29|) is nothing but the total angular momentum, L, 
of the system (with respect to the center of mass): 

N 

Ju^av = E ^fc^fc = L. (34) 
k=l 

It is thus analogous to the formula 

N 

mR = E '^k^k = P (35) 

k=l 

for the total linear momentum, p, of the system. Using the linear and angular 
momenta the total kinetic energy of a system can now be written 

Here m is the total mass, which is always constant, and p the total linear 
momentum, which is constant in the absence of external force. Thus in the 
absence of a net external force on the system the first term in this expression 
for T is constant. What about the second term? In the absence of a net 
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external moment (of force) on the system L is constant, but the (inverse) 
inertia tensor is, in general, not since it depends on the 'shape' of the 
system. 

In an isolated body, such as a star or planet, one can expect that internal 
dissipative forces, in the long run, will make the internal relative motions zero 
so that T" — > 0. Such bodies will thus end up having only center of mass 
translational and average rotational energy and both will then be constant. 

8 Conclusions 

The fixed axis formalism at the beginning of this paper is simple and useful 
enough to be included in even fairly elementary texts, though I believe it 
is presented here for the first time in the international literature. I have 
introduced part of it in course notes that are used in the school of engineering 
physics at KTH. The point of view that the Jellinek and Li (1989) separation 
of a rotational part of the kinetic energy is analogous to the well known 
separation of center of mass translational energy, is pedagogically useful at 
higher levels, and the general ideas deserve to be better known. It is the 
opinion of the author that most of the equations and results of this paper, 
in fact, belong in a comprehensive advanced mechanics course. 
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